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Abstract 

We study the physics potential of a neutrino superbeam experiment with a 2540 km baseline. We 
assume a neutrino beam similar to the NuMI beam in medium energy configuration. We consider 
a 100 kton totally active scintillator detector at a 7 mr off-axis location. We find that such a 
configuration has outstanding hierarchy discriminating capability. In conjunction with the data 
from the present reactor neutrino experiments, it can determine the neutrino mass hierarchy at 3a 
level in less than 5 years, if sin 2 2#i3 > 0.01, running in the neutrino mode alone. As a stand alone 
experiment, with a 5 year neutrino run and a 5 year anti- neutrino run, it can determine non-zero 
6»i3 at 3a level if sin 2 26 13 > 7 x 10~ 3 and hierarchy at 3a level if sin 2 20 i3 > 8 x 10~ 3 . This data 
can also distinguish Sep = tt/2 from the CP conserving values of and 7r, for sin 2 29n > 0.02. 
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I. INTRODUCTION 



During the past decade, great progress has been made in the study of neutrinos. In 
particular, the mass-squared differences and two of the mixing angles have been measured 
with good precision. Despite this progress, the following properties of neutrinos are still 
unknown 

• The absolute mass scale of the neutrinos, 

• Dirac vs Majorana nature of neutrinos, 

• The value of the mixing angle #13, 

• Mass pattern or mass hierarchy of the neutrinos, 

• CP violation in neutrino sector and the value of the CP- violating phase Sep, 

• Deviation of #23 from 7r/4. 

KATRIN tritium decay experiment and various neutrinoless double beta decay 

experiments [3j currently underway will address the first two problems. The next four 
problems fall under the purview of neutrino oscillations and hence are expected to be solved 
by the current and future neutrino oscillation experiments. There are three efforts underway 
to determine the mixing angle #13 using reactor neutrino data. They are (a) Double Chooz, 
expected to start taking data this year Q with sensitivity In' 2ft, > 0.04, (b) RENO fl 
with sensitivity sin 2 2# 13 > 0.02, and (c) Daya Bay |6[ with sensitivity sin 2 2# 13 > 0.01. 
In addition, there are two accelerator experiments, T2K and NOz/A, which can determine 

ri 

non-zero #13. Their sensitivities are sin 2#is > 0.01 at 90% confidence level for T2K m and 



; 13- 

sin 2 26 13 > 0.03 at 3a for NOz/A 



,8 



Thus it is hoped that the determination of #13 can be 



achieved by these experiments within the next decade. 

However, it was shown recently 9] that the currently running or planned experiments will 
be able to determine the mass hierarchy or Sep only for very favourable values of #13 and Sep- 



The effect of mass hierarchy and Sep are closely interlinked 10Ml4l|. Disentangling these 



m 



16| 



two effects requires considerable effort and possibly data from many experiments 
Thus the question of designing a neutrino beam-detector configuration, which will enable us 
to determine hierarchy and Sep unambiguously, acquires great importance. 
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There have been numerous proposals made to achieve the above purpose. These include 



megaton size detectors 



17H19|| , a pair of detectors with different baselines 



3ut with the same 



value of L/E 15|, |l6| and detectors at the 'magic' baseline of 7500 km 



20 



2lj. Recently 



we made a proposal of a shorter magical baseline [22J, which has exceptional capability to 
determine mass hierarchy, independently of Sep- In this report, we study the full physics 
potential of our proposal. 



II. CALCULATION 



As mentioned before, the unknowns to be determined from the neutrino oscillation ex- 
periments are the mixing angle #13, the neutrino mass hierarchy, the CP violating phase 
5c p and the octant of 6 23 . The — )■ v e oscillation probability P Me is sensitive to all these 
unknowns. In this paper, we will not consider the problem of octant ambiguity of 6 23 . We 
will assume the true value of 6*23 to be its best fit value t t 1 4 . The expression for P^ e , in three 
flavour oscillations including matter effects, is given as 

sin 2 ((l - A) A) 
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where A = (1.27A 31 L/E), A = A/A 31 and a = A 21 /A 31 w ±0.04, as given by the solar and 
atmospheric neutrino data [26]. Here L is in km, E is in GeV and A31 is given in units of 
eV 2 . The coefficients, Ci are given by 



C = sin 2 #23 sin 2 26 13 , 

C\ = cos 6*i3 sin 2^12 8^26*13 sin 2^3 cos(A + 5cp) 
C 2 = sin 2 26 12 cos 2 6 23 . 



(2) 
(3) 
(4) 

The matter effects are given by the matter term [27j A (in eV 2 ) = 0.76 x 
10~ 4 p (gm/cc) E (GeV) eV 2 . p is the density of the matter through which the neu- 
trino propagates. A31 and hence A are positive for normal hierarchy (NH) and negative 
for inverted hierarchy (IH), for both neutrinos and anti-neutrinos. A, on the other hand, is 
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positive for neutrinos and negative for anti-neutrinos. Therefore, A is positive for NH and 
negative for IH in the case of neutrinos, and vice versa for anti-neutrinos. 

The dependence of P Me on Sep and 9 13 is seen explicitly in Eq. (pQ). In addition, the 
quantities A and A are sensitive to the neutrino mass hierarchy. Thus, this channel which 
is sensitive to all the three unknowns mentioned above, is a good candidate to probe these 
unknown parameters. From Eq. ([T]) we see that P Me is large for A = (1.27 A^L/E) ~ 7r/2. 
On the other hand, the matter term A becomes significant at large energies. Therefore, the 
matter term starts to play an important role in shaping the oscillation probabilities only for 
long baselines. For a baseline of 732 km, and for fixed values of #13 and Sep, it increases P^c 
by about 25% for NH and decreases it by the same amount for IH 28) . For longer baselines, 
the change in the probability is even larger. At shorter baselines, the sensitivity of P^e to 
the matter term and hence the mass hierarchy, is low. 

Since P Me depends on three unknowns, 613, mass hierarchy and Sep, it is impossible to 
determine any one of them by making a single measurement. Evenjf one measures P Me 



as a function of energy, it is possible to get degenerate solutions [10Hl4|. For example, 
it is possible to have P^ e (6 , i 3 , NH; E) = P^ c (9[ 3 , IH; E), as illustrated in Fig [TJ Such a 
degeneracy can be broken, if the value of #13 can be determined from different experiments. 
Reactor anti-neutrino experiments with L 1 km, can determine #13 unambiguously. The 
currently planned experiments can measure non-zero values as small as sin 2 2# 13 > 0.01 
4-6]. Thus we expect that this degeneracy can be lifted by the results of reactor neutrino 
experiments. 

Even if the ambiguity due to the value of #13 is lifted, there is a further ambiguity between 
hierarchy and 5c p- This ambiguity arises due to the following reason. If we replace v — > T>, 
we have 5 c p — > —Sep and A — » — A. However, if we flip the hierarchy, we have A — > — A 
and A — > —A. Hence, the change in P^ e due to hierarchy flip can be compensated by 
changing Sep- This ambiguity is particularly important for sin 2 2#i3 < 0.02. Therefore, 
even if we know the value of 9\s, it is possible to get degenerate solutions of the kind 
P^ C (NH, 5cp; E) = P^ e (IH, 5' CP ; E) over a wide range of energies, as illustrated in Fig. |2j 
We are led to consider various improvisations to overcome this problem. 

In Eq. (TjQ), only the C\ term contains 5c p- The hierarchy-<5cp ambiguity can be lifted 
if the term containing C\ can be chosen to be zero. One simple way to do it is to choose 
A A = 7r [29] . This condition can be achieved for both NH and IH simultaneously and gives us 
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FIG. 1: Hierarchy-6>i3 degeneracy in P MO at 1300km. We see that P, IC ( NH, sin 2 29 13 = 0.02; E) = 
P^ C {IH, sin 2 26»i3 = 0.05; E). The value of S C p is taken to be 330° here. 
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FIG. 2: Hierarchy-6 C p degeneracy in P MC at 1300km. We see that P^NH, S C p = 45°; E) = P^IH, S C p = 
225°; E) for £ > 1.5 GeV. The value of sin 2 29 13 is taken to be 0.02 here. 



an energy independent condition L m 7500 km. This is the famous magic baseline condition 
and it has been extensively studied 20, 12I]. In 22] it was noted that the C\ term can be 
made zero for IH only by choosing (1 — A) A = —71. This makes P jUe (Ji?) independent of Sep 
and also very small (< a 2 m 0.002). If we choose (l-i)A = vr/2 for NH, then P lie (NH) will 
be quite siginificant for all values of Sep- Solving the above pair of simultaneous equations, 
we obtain L = 2540 km and E = 3.3 GeV. For this baseline, at this energy, there is a clear 
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separation between P^ e (NH) and P^ C (IH) as illustrated in Fig. [3J This is true for values 
of #13 at least as small as those detectable by the current reactor experiments. In addition, 
P^ci^NH) retains its dependence on Sep unlike in the magic baseline case. This, in turn, 



can be helpful in measuring 5c p 



30]. 



In our calculations, we have considered a superbeam experiment with a baseline of 2540 



18| and 



km. Note that the distance from Brookhaven Laboratory to Homestake mine [l7, 
that from CERN to Pyhasalmi 31] are close to our shorter 'magical' baseline. Since our 
energy range of interest is around 3 — 4 GeV, the neutrino source must have unoscillated 
event spectrum peaking in this energy range. The NuMI beam in the medium energy option 
has an unoscillated event spectrum that peaks at 3.5 GeV for locations 7 mr off the beam 
axis js]. 



III. RESULTS 



We have done all our calculations using the software package GLoBES 32|, |33[. The 
source is taken to be a NuMI like beam in the medium energy configuration. The detector 
is assumed to be a 100 kton totally active scintillator detector whose capabilities are similar 
to NOz/A. Such a detector is placed at a distance 2540 km from the mr off-axis 

location. The fluxes for this location were calculated using the program from [34]. The beam 
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power is taken to be such that it corresponds to 10 x 10 20 POT/yr. This is about 30% more 
than the design power for NOz/A. Thus we have an exposure of 1000 x 10 20 kton-POT/yr. 
This is about 10 times larger than the exposure for NOz/A. The baseline in our proposal is 
thrice as long as that of NOz/A. Hence the factor 10 increase in exposure makes the data to 
have the same statistical weight as that of NOz/A. 

We have computed the v e appearance spectrum in the detector for the energy range from 
250 MeV to 10 GeV, in bins of 250 MeV width. A Gaussian energy smearing function 



with a width of O.lyE (GeV) has been assumed for these events. We have considered four 
different types of events in our analysis: signal (z/^ — > z/ e ) events, beam background events 
(that is, events caused by intrinsic v e component of the beam), background events due to 
misidentified muons and the background events due to neutral current interactions. Energy 
independent cuts identical to those for NOz/A have been imposed on these channels, which 
substantially enhance the signal to background ratio 0, [j^ . 



In our calculations, the so 



ar parameters sin 2 9u = 0.304 and A21 = 7.65 x 10 5 eV 2 have 



been kept fixed throughout |26|, ]36|, [37J. The values of A 3 i and 6*23 are taken from three 



flavour oscillation fit of the atmospheric 38] and MINOS data 391 ]. We take the true value 
of sin 2 2#23 to be its best fit value 1. We assumed an error of 2% on this measurement. Hence 
the test values of #23 are to be chosen to satisfy the constraint 0.94 < sin 2 2# 2 3 < 1- Thus we 
have the range 38° < #23 < 52°. The situation regarding A 31 is a bit more complicated. The 
most accurate measurement of the larger mass-squared difference is given by the MINOS 
experiment 39] . If the results of MINOS are interpreted in terms of three flavour oscillations, 
the effective mass-squared difference measured by it is not |A 3 i| but is the magnitude of a 
linear combination of A 31 and A 2 i 40, 41]. It is shown j^l ] that the linear combination is 



of the form A 31 + A 2 if{9 i j ) 5cp), where 



5cp) = cos 2 612 - cos 5 C p sin6 13 sin26>i 2 tan6> 2 3- 



(5) 



In computing P Aie (NH) and P Me (IH) we must relate the corresponding values of A 31 via 

A 3 i (IH) = -A31 (NH) - 2A 2 i/(%, 5cp), (6) 

rather than as A 3i (IH) = -A 3l (NH). If we were to take A 3i (IH) = -A 3l (NH), then we 
would get a significant fake hierarchy sensitivity in the z/ M disappearance channel, even for 
very small values of 6*13. In our calculations, we have related A 31 (NH) and A 31 (IH) via 



Eq.®, and used the best fit value A 3i (NH) = 2.4 x lO~ 3 eV 2 , with an error of 5%. The 
systematic errors include a 5% normalization error and 2.5% tilt error on both the signal 
and the background. We have also taken into consideration a 5% error in the density profile 
of the earth. 

First we computed the hierarchy sensitivity of this setup, running it in the neutrino mode 
only. In this part of the calculation, we limited ourselves to the range 0.01 < sin 2 2#i3 < 0.1. 
That is, we assumed that this parameter will be determined by one of the current 
experiments. First we assumed NH to be the true hierarchy and computed the x 2 f° r 
distinguishing this from IH, as a function of sin 2 26 ) i3. In computing the x 2 ; we combined the 
results of our proposed set up with those of the reactor experiments and of T2K and NOz/A 
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4J4|. The data from the reactor experiments provide a 'prior' on the values of 
sin 2 26 l i3. In computing x 2 , we have marginalized over the range 0.01 < sin 2 26*13 < 0.1. We 
have varied the CP phase over the full range — 71 < 5c p < tt in computing the event rates 
both for NH and for IH. We have then repeated the above calculation, this time assuming 
IH to be the true hierarchy. The results of these calculations are tabulated in Table [H We 
indicate here the exposure time (in years) required for hierarchy discrimination up to a 3<r 

13- 



confidence level function of sin 28 



From the results of the table, it is worth remarking that, the present 15 kton NOz/A 
detector, with its current design luminosity, can determine the neutrino mass hierarchy, 
independently of 5cp, in a five year neutrino run, if sin 2 26*13 > 0.03, were the detector 
placed at a distance of 2540 km with an appropriately designed beamline. 

Next, we have calculated the ability of our setup (in conjunction with T2K, NOz/A and the 
reactor experiments) to determine 5c p as a function of sin 2 26*13. in addition to the data from 
a five year neutrino run, we assume data from a five year anti-neutrino run. We have chosen 
fifteen different pairs of true values of sin 26*13 and 5c p given by sin 2 2^13 = 0.01,0.05,0.10 
and 5cp = —160°, —80°, 0°, 80°, 160°. We assumed normal hierarchy and computed the 
allowed la and 2a regions in the sin 2 26*13 — 5cp plane. In Figure HI we have plotted the 
outer contours of these allowed regions, corresponding to the fifteen different true parameter 
sets. The panel on the left gives the results including the systematic errors, whereas the 
panel on the right shows the allowed contours without the systematic errors. From the right 
panel, we see that the allowed regions are very close to the true values and the value of 
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sin 2 26> 13 (true) 


Exposure time(NH) 


Exposure time(IH) 


0.10 


0.022 


0.048 


0.09 


0.026 


0.057 


0.08 


0.031 


0.068 


0.07 


0.040 


0.082 


0.06 


0.051 


0.105 


0.05 


0.070 


0.137 


0.04 


0.104 


0.195 


0.03 


0.180 


0.420 


0.02 


0.425 


2.600 


0.01 


2.950 


4.800 



TABLE I: Exposure time in years required for 3a hierarchy discrimination, for 1000 x 10 20 kton- 
POT/yr. The second (third) column shows the results if NH (IH) is the true hierarchy. It is 
assumed that the experiment runs in neutrino mode only. 

5cp can be determined with an uncertainty of about 10°. The value of sin 2 2#i 3 is rather 
tightly constrained by the data from reactor neutrino experiments. However, the inclusion 
of systematic errors makes the allowed regions considerably larger, sin 2 2# 13 is allowed to 
deviate from its true value by about 0.005 and the uncertainty in determining Sep rises to 
about 20°. 

We have also computed the physics capabilities of the proposed setup as a stand alone 
experiment. As mentioned before, in the previous calculations, we have taken a lower limit 
of 0.01 < sin 2 2#i3. The reason for this is, that for values of #13 above this limit, the data 
from the reactor experiments provide a 'prior'. In the calculations discussed below, we did 
not impose this lower limit on #13 but allowed it to be as small as sin 2 29 13 ~ 10~ 4 . 

In computing the results discussed below, we have assumed a 5 year neutrino and a 5 
year anti-neutrino run. Based on this data, we computed the sensitivity of this set up to 

• measure non-zero 9 13 

• determine hierarchy 

9 
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FIG. 4: la and 2a contours for sensitivity in the sin 2 26>i3 — Sep plane. The plot on the right is without 
systematics. 

• rule out the CP conserving value of Sep = or tt. 

In figure [5], we plot the true values of sin 2 2# 13 for which this parameter can be distin- 
guished from 0. We do this as a function of the true value of Sep for both NH and IH. If 
NH is the true hierarchy, then values of sin 2 2#i3 as small as 8 x 10~ 4 can be measured at 3a 
level, for advantageous values of 5cp (between — n/3 to — 2tt/3). For adverse values of 5cp 
(between n/6 to vr/2), however, #13 can be distinguished from only for sin 2 29^ > 7 x 10~ 3 . 
If IH is the true hierarchy, the results are similar, except for the fact that the advantageous 
and adverse values of 5c p are exchanged. 

In figure [6j we plot the sensitivity of this setup for distinguishing the mass hierarchy as a 
function of the true value of sin 2 2#i3 and Sep- Once again, for advantageous values of Sep, 
the mass hierarchy can be determined at 3a for sin 2 2#i3 > 8 x 10~ 4 , whereas for adverse 
values, it can be determined only for values which are about an order of magnitude larger. 

Here we must insert a note of caution. Figures [5] and [6] give the sensitivity of our proposed 
setup as a function of the true value of Sep- However, the true value of Sep is unlikely to 
be known from any of the current experiments. Hence, while making a statement regarding 
the sensitivity of our setup, we must consider the most advserse value of Sep- Thus, a five 
year neutrino and a five year anti-neutrino run in our setup can measure, independently of 
Sep, sin 2 2#i3 if it is larger than 7 x 10~ 3 and it can determine hierarchy for a similar lower 
limit on 13 . 

Finally, figure [7] shows the ability of the setup to exclude the CP-conserving case, Sep = 

10 



0, 7T. We find that the maximal phase Sep = ±7r/2 can be distinguished from 5cp = 0, 7r for 
sin 2 2#i3 > 0.01. But for other values, discrimination from CP conserving case is possible 
only for larger values of #13. In fact, for 5cp outside the two ranges — 2ir/3 to — 7r/3 and 
7r/3 to 27r/3, the discrimination from CP conserving case is not possible even for the largest 
allowed value of #13. 

H2B collaboration 45) have considered a neutrino superbeam experiment with a baseline 



in excess of 2000 km. Their setup envisaged directing a 4 MW neutrino superbeam, with peak 
flux at E v = 4 GeV, from Tokai in Japan to a 100 kton water Cerenkov detector at Beijing, 
2100 km away. Such a set up is close to satisfiying the shorter magic baseline conditions 



derived in section-II. They found that [46|, |47) such a setup, running only in neutrino mode 
for five years, can determine mass hierarchy for all values of Sep, if sin 2 26* 13 > 0.02. 
There was a previous study of a superbeam proposal with 2540 km baseline 
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18). 



In this proposal, a wide band neutrino beam, with 0.5 < E u < 7 GeV energy range from 
Brookhaven National Laboratory is directed at a 500 kton water Cerenkov detector at Home- 
stake mine. The beam is expected to be produced by 28 GeV protons with 1 MW intensity. 
This corresponds to 22 x 10 20 POT/yr assuming one year contains 10 sec running time 



49) . The detector is assumed to be along the beamline, rather than at an off-axis lo- 



cation, which leads to the rather wide range in the energy. With a 5 year run in neutrino 
mode only, this setup can set a 90% confidence level upper limit sin 2 26* 13 < 0.005 if NH is 
the true hierarchy and sin 2 26*13 < 0.02 if IH is the true hierarchy. It can also determine 5cp 
with about a 25% uncertainty. 

A few years ago, a US long baseline neutrino study group advocated a 1300 km baseline 
experiment to determine the unknown neutrino parameters from oscillations 48, 4jJ. They 



considered a wide band neutrino beam with 0.25 < E u < 7 GeV energy range directed 
from Fermilab to a 300 kton water Cerenkov detector at Homestake mine. They assumed 
a neutrino run with 30 x 10 20 POT with a similar anti-neutrino run. The exposure of this 
setup, in POT-kton, is twice that of the setup we considered. The CP-violation discovery 
potential of this setup is better but the hierarchy determination potential of our proposal is 
better Q- 

Recently, there was another study which considered additional magical properties of the 
2540 km baseline [3^]. In that study, the following set up was considered. The neutrino 
source is a low energy neutrino factory with a flux of 5 x 10 21 muons/year with the muon 
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energy of 5 GeV. The detector is a 25 kton totally active scintillator detector and it is 
assumed that this setup will run with positive muons for 2.5 years. With these beam and 
detector specifications, the authors showed that #13 can be measured for values sin 2 29 13 > 
4 x 10 -3 and the hierarchy can be determined for sin 2 2# 13 > 8 x 10~ 3 . This set up seems to 
have a good ability to discover non-zero Sep- 

180 
120 
60 

-60 
-120 

-180 

1 

FIG. 5: 3cr and 5<r contours in the sin 2 2#i3 — Sep plane for excluding #13 = 0. 



IV. CONCLUSION 

In an earlier paper, we demonstrated that a neutrino superbeam experiment with a 2540 
km baseline is particularly well suited to determine neutrino mass hierarchy, independently 
of the value of 5c p- In this paper, we performed a full-fledged study of the physics potential 
of such an experiment. We considered a 100 kton totally active scintillator (NOz/A like) 
detector with the neutrino source being a NuMI-like beam in medium energy option with 
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FIG. 6: 3cr and 5a contours in the sin 2 29r3 — Sep plane for excluding the 'wrong' hierarchy. 




FIG. 7: 3<7 and 5a contours in the sin 2 2#i3 — Sep plane for excluding CP conservation. 

10 x 10 20 POT/yr. We performed the analysis using GLoBES software, where all the appro- 
priate background events have been included and background suppression factors relevant 
for NOz/A are imposed. We demonstrate that a five year neutrino run in this setup, together 
with data from reactor neutrino experiments, can determine the neutrino hierarchy inde- 
pendently of Sep for sin 2 26 ri > 0.01. This statement holds irrespective of whether the true 
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hierarchy is NH or IH. With additional data from a five year anti-neutrino run, this setup 
is capable of determining 5c p with an uncertainty of about 20°. 

If this setup is considered by itself, then the data from a five year neutrino run plus a 
five year anti-neutrino run can measure non-zero #13 and determine neutrino mass hierarchy 
for sin 2 2# 13 > 7 x 10~ 3 . These statements hold true for all possible values of Sep- It is also 
possible to obtain signals for CP violation, for moderately large values of sin 2 2^13. 
Acknowledgement We thank Ravi Shanker Singh for collaboration on an earlier paper 
where the magical properties of 2540 km baseline were first discussed. We thank Srubabati 
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